International Journal of Advanced Engineering Research and Science (IJAERS) 
https: //dx.doi.ora/10.22161/iiaers.5.12.25 


[Vol-S, Issue-12, Dec- 2018] 
ISSN: 2349-6495(P) / 2456-1908(0) 


Fractional Order Butterworth Filter 

Mehmet Emir Koksal 


Department of Mathematics, Ondokuz Mayis University, 55139 Samsun, Turkey 
Email: mekoksal@omu.edu.tr 


Abstract —Fractional order Butterworth filter derived 
from the conventional 3 rd order Butterworth filter by the 
transition of the ordinary derivative to the fractional one 
as in [1] is investigated in this paper. The change of the 
filter characteristics is studied depending on the order of 
the fractional derivative. Effects of the transformation on 
the components’ behaviors of the filter is formulated. 
Design of the filter in the sense of choosing the filter 
parameters to satisfying the given specifications is 
described. 

Keywords—Butterworth filter, control system, fractional 
order. 


I. INTRODUCTION 

Fractional calculus (FC) [2,3]provides more accurate 
models than classical calculus for many of the physical 
systems showing intrinsic fractional order (FO) 
behavior[4-6]. Fractional derivative takes into account the 
past history when memory like behaviors exhibit in 
components [7-10]. The same advantages appear in some 
electrical circuits as well [11-16]. 

In [1], a method is given for obtaining fractional order 
differential equation (FODE) from ordinary differential 
equation (ODE) keeping the consistency of 
dimensionality in physical systems. This method is used 
to find the natural response of a FO RLC circuit by F. 
Gomez at al. [17]. 

In this paper, FO Butterworth filter (FOBF) is derived by 
using a method similar to Gomez’s method. The ordinary 
derivative is replaced by FO derivative) by the 

transformation 

d 1 d Y 

- * —-, ( 1 ) 

dt a 1 7 dt Y 

wherey is an arbitrary parameter defining the order of 

derivative and the auxiliary parameter a has the 

dimension of second. Time and frequency response 

characteristics of the obtained FOBF is investigated 
according to the transformation parameters yand er. 

The paper is organized as follows: hi Section 2, a brief 
summary of the FC relevant to the content of this paper is 
presented. Section 3 covers the short introduction of the 
conventional Butterworth filter. Section 4 describes the 
evaluation of the FOBF and the related formulas. In 
Section 5, the time and frequency domain characteristics 


of the FOBF is investigated. Finally, conclusions are 
included in Section 6. 


II. FRACTIONAL CALCULUS 


It is well known that Caputo fractional time derivative 
[16, 17] of a function fit) is defined by 

- W'ta 

r(n — y)J 0 (t — T) Y ~ n+1 V J 
wherer EJ?,n-l<y<nelV = {l,2,...}, f n 

represents the ordinary derivative of order n , and T is the 
Gamma function. For y = 1, Eq. (2) gives the usual 
derivative. It is assumed in the scope of this paper that t 
represents the time in seconds (s), and n = 1 so that (2) 
reduces to 


d Y 

dtyf 1 ~ T(1 

We note that Eq. (2b) yields 
d Y rs 

- fit ) 

dt yJ 


- r)Jo 


- V) Jn 


rco 

it - r) Y 

0 f n CO 


dr. 


(2b) 


dr = 0. (2c) 


t=o r(i-y)J 0 (t-ry 

Laplace transform of Eq. (2a) is 
£{ C o D Y t fit)} 

71-1 

= s y Fis ) — ^ s Y ~ k ~ 1 f ^ (0 + ), n — 1 < y < n. (3a) 


For n = 1. that is 0 < y < 1, 

L { C qD y fit)} = s Y Fis) — s Y ~ 1 fiO + ), 0 < y < 1. (3b) 
Considering the inverse Laplace transform, a few useful 
formulas are fisted as follows: 

L _1 = ( 4a ) 

Where £ a pi,) is the 2-parameter generalization of the 
Mittag-Leffler function. It is defined by 

00 k 

=Z r iB + akr CC ' P ’ Z£C)Re< ' a) > °' 

Re(p) > 0, (4b) 

which reduces to the original Mittag-Leffler function for 
0=1: 

£ a iz ) = Sa.liz) 

00 k 

= Erd + akV 
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III. LOW PASS BUTTER WORTH FILTER 

The transfer function of the 3 rd order low pass 
Butterworth filter (LPBF) normalized with the 3 dB 
cutoff frequency of a> 0 =1 r/s is given by 

(s + l)(s 2 + s + l) 

1 

(s 3 + 2s 2 + 2s + 1) 

One of the circuit realizations of this transfer function is 
done as the voltage ratio transfer function as it is shown 
in Fig.l. The analysis of the circuit results the following 
transfer function: 

Ms) 


Vs (s) 


His) 


o R r j L- 1 +Lo R 

s 3 + —s 2 + ——~ s +- 

Z<2 l^2 c L ± L 2 C 


(5b) 


Equating the coefficients of the transfer functions in Eqs. 
(5a) and (5b) and using the normalized resistance R = 
1 fl, the following component values are obtained for the 
filter. 

4 

L 1 = 1.5 H, L 2 = 0.5 H, C = - F, R = 1 fl. (5c) 



+ 

v rW 


Fig.l: 3 rd order Butterworth filter with o) 0 = 1 r/s. 

The differential equation relating the input voltage v s (t) 
to the output voltage v R (t) is obtained formEq. (5c) by 
inserting the Laplace operator s with the derivative 
operator d/dt and arranging the temis; the result is 
d 3 „ , R d 2 , . T 1 +T 2 d 


;v R (t) + 
R 

+ - 


L 2 dt 2 


v R (t) + 


L t L 2 C 


v R (t) 


R 


L 1 L 2 C 


L t L 2 C dt 

MO- 


? R (t) 


( 6 ) 


The frequency domain characteristics, namely gain and 
phase responses of the filter are obtained replacing 
s by ja> in Eq. (5a); the result is 

H(jax) = M(a>) eM "' 1 ,where (7a) 

M(m) = I Hi] ox) I = , 1 , (7b) 

Vl + O ) 6 

. . _ (2 - co 2 )co 

0(a>) = ArgiHijo))} =-tan 1 - — ^ 2 , (7c) 

are the phase and phase characteristic, respectively. From 
Eq. (7a), it is obvious that the gain decreases from its 
value 1 at a> = 0, and it becomes equal to its half power 
value M = 1 /s[2 at cu = 1, it decays to zero as a> —> 00 . 
And the phase start from 0 at a> = 0, it decreases 
monotonically to —270° as a> -» 00 . The time and 
frequency characteristics are not plotted at this stage since 


they will appear as the special case of the fractional 
Butterworth filter considered in the next section. 

IV. FRACTIONAL ORDER LOW PASS 
BUTTERWORTH FILTER 

Applying the fractional transformation to the components 

of the LPBF in Fig. 1, we have the following time and 

Laplace domain element behavior equations in the circuit: 

v R (t) =Ri R (t), V R (s) = RI r (s), (8a) 

. . 1 d r i L (t) 

v,.(t) = L:— - - -,17,. (s) 

Ll l a 1 ~ y dtr Ll 

1 

= L — —sVI L .(s), i = 1,2, (8b) 


ir(t) 


l o 1 -r 
1 d y v r (t) 


1 = c —filZ^Vc (s) ■ (8c) 


c “ a dt y 
The frequency and time domain analysis of the circuit 
with these component behavior equations, or direct 

in Eq. (5b) and — —> 


substitutions s 

1 d y d 2 


(T 1 ~y 
1 d 2y d 3 


d 3y 


a'-Y dt y ' dt 2 a 2 0 --Y) dt 2 Y ’ dt 3 a 3 H~y) dt 3 Y 

yield 

V R (s) 


■, in Eq. (6) 
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dt 3y K L 2 dt 2y 

(L 1 +L 2 )cr 2(l_r ' ) d y f ^ 
+ . . _ - ~v R (t) 


L 1 L 2 C 
RcjIH-y) 


+ ■ 


L 1 L 2 C 


v R (t) = 


dt y 

Ra^-Y) 


L 1 L 2 C 


-v s (t). (9b) 


respectively. The characteristic equation is obtained by 
equating the denominator polynomial of the transfer 
function in Eq. (9a) to zero: 

Rc t ( M M+MMM Ra 3(l - y) 


s 3y + 


; 2y 


+ ■ 


■S^ + - 


L 2 L 1 L 2 C L 1 L 2 C 

= 0. (9c) 

The characteristic polynomial appearing in the left side of 
the equality in Eq. (9c) is s a commensurate polynomial 
in power s y . So the time domain responses (such as step 
and impulse responses) of the filter can be obtained 
analytically by using Mittag-Leffler function [5]. But, 
these solutions are not included within the content of the 
paper and it is satisfied by their plots only, instead we are 
confined to the frequency domain responses. 

V. FILTER CHARACTERISTICS 

To find the gain and phase characteristics of the FO LPBF 
derived in the previous section, we let s = ja> in Eq. (9a) 
and use the identity 
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Y p >2 T 


a> r e 2 


r TL\ y 

(ja>y = 0>y [e ] 2) = 

= ^[cos^y) + jsin (—y)] 


The result is 


M(a>) = 


0(u>) 




Vi4 2 + B 2 ’ 
(B\ 

Arctan — , 

\aJ 


where 


u> 3y cos(l.5ny) + a 2 a> 2y cos(ny ) 


2 a> 2y cos (ny) 


+a 1 (i) r cos(0.5ny) + a 0 , 

B = a> 3y sin(l.57ry) + a- 

+a 1 a> Y cos(0.57ry), 
R(j^- y) (L 1 +L 2 )o- 2(1_r) 


L 1 L 2 C 


Ro 3 ^ 


L 1 L 2 C 


( 10 ) 

(Ha) 

(lib) 
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(lid) 
(He, f) 
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The gain and phase characteristics on the logarithmic 
scale (Bode plots) of the filter for a = 1 and different 
values of y is shown in Fig. 2 where the conventional 3 rd 
order Butterworth filter characteristics is shown by the 


thick dashed red line (y = 1); The details of the 
numerical data is given in Table 1 where M p (in dB) is the 
peak gain at the peak frequency a> p (in rad/s) , M 1 = 
M 2 = M p — 3are the gains at the cut off frequencies a> 1 
and o) 2 , BW is the bandwidth which is defined as ui 2 — 
uq if a peak exist (the case y = 1); otherwise BW = a> 2 
(the cases y = 1, 0.75, 0.50) The quality factor Q = 
m 2 /(m 2 — &q)is defined only for y = 1.25 for which the 
gain characteristics has a peak exceeding the 0 dB level 
and the filter can be considered as of bandpass type as 
well. It is observed from Fig. 2 and the data given in 
Table 1 that the first cut off frequency &q and the 
associated gain M 1? and Q are defined only when the 
characteristic is handled as a BP type for which 
BW=0.221. In general, the bandwidth (BW) and hence 
the cut off frequency (u> 2 = BW) decreases with 
decreasing y values for the LP filter. 

The phase characteristics decrease from 0° to — 180° as 
a>: 0 -» oo but with a faster rate at the intermediate 
frequencies with increasing y. 
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Bode Diagram 
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Fig.2: Bode plots of the transfer fimction for a =1.00. 
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Table.l: Numerical data for frequency response curves when a = 1 for the filter. 


Y 

mi - (MJ 

M 2 ) 

- (M P ) 

BW - (BP) 

Q - BP 

1.25 

0.852 - 8.841 

1.073 - 8.867 

0.979 -(11.844) 

1.073 - (0.221) 

4.430 

1.00 

— 

1.000 - (-3) 

0- (0) 

1.000 - (—) 

— 

0.75 

— 

0.313 - (-3) 

0- (0) 

0.313 - (—) 

— 

0.50 

— 

0.058 - (-3) 

0- (0) 

0.058 - (—) 

— 


The step response of the filter is shown in Fig. 3 for 
values of = 1.25,1.00,0.75,0.50 . It is seen that all the 
responses start from 0 at t = 0 and approaches to 
oo as t -» oo, which is a typical property for a LP filter. 


The response is highly oscillatory for y = 1.25, the case 
in which the filter can be interpreted as a BP filter as well. 
There is a smaller overshoot for the ordinary 3 rd order 
Butterworth filter (case y = 1.00, thick red dashed line). 
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The rise time and the settling time decreases as y takes 
smaller values. See Table 2 for the detailed numerical 
data, where T r : rise time, T r : rise time, T pl , T p2 , T p3 : 
peak times, M pl ,M , M p3 \ peak values, T s : settling time all 


in seconds. Some of the data could not be detected since 
the simulations are done up to 40 s. Rise time increases 
with decreasing/ and the minimum settling time occurs 
for the conventional filter (y). 



Fig. 3: Step responses of the transfer junction when a =1.00. 


Table.2: Step response characteristics for a =1.00. 


Y 

T 

i r 

Tpi M pl 

Tp2 M p2 

Tp 3 Mp3 

T 

1 s 

1.25 

1.49 

4.56 - 1.678 

10.95 - 1.315 

17.27 - 1.634 

24.63 

1.00 

2.29 

4.93 - 1.082 

12.10 - 1.002 


5.98 

0.75 

11.25 




27.51 

0.50 







Similar characteristics for a = 0.1 and the same different 
values of y considered before are presented in Fig. 4 and 
Table 3. The general arguments discussed for the case 


affected. For y < 1. all the critical frequencies decrease 
with decreasing y. 

Considering the phase characteristic, it is slightly 


cr = l hold. For y = 1.25, all the critical frequencies are 
increased whilst the critical gains and the quality factor 
remain the same. For y = 1.00 the responses are hardly 


increased (decreased) for y = 1.25 (y < 1.00) and almost 
unaffected for y = 1.00. 


Bode Diagram 
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Table.3: Numerical data for frequency response curx’es when o = 0.1 for the filter. 


Y 

mi - (AfJ 

W 2 iM \2 ) 

a> p - (. M p ) 

BW - ( BP ) 

Q - BP 

1.25 

1.350 - 8.835 

1.701 - 8.854 

1.552 -(11.844) 

1.701 -(0.351) 

4.422 

1.00 

— 

1.000 - (-3) 

0- (0) 

1.000 - (—) 

— 

0.75 

— 

0.145 -(-3) 

0- (0) 

0.145 -(—) 

— 

0.50 

— 

0.006 - (-3) 

0-CO) 

0.006 - (--) 

— 


Step response of the filter for er = 0.1 is shown in Fig. 5. 
It is seen that the conventional Butterworth filter response 
is not affected by a when y = 1.00 (see the legends y = 
1.25 and y = 1.25 *). Further, the peak values of 
oscillations are not affected by changing erfrom 1 to 0.1. 


But there is an apparent time lead (squeeze) for y = 
1.25 >, and time lag (spread) for y = 0.75 < 1. The lag 
(spread) is higher for y = 0.50. See Table 4 for numerical 
details. The time lead (lag) is used in the sense of faster 
(slower) motion. 



Fig. 5: Step responses of the transfer Junction when er =0.1. 


Table 4: Step response characteristics for er =0.1. 


Y 

T 

i r 

Tpi Mpi 

Tp2 M p2 

7p>3 M p3 

T 

1 s 

1.25 

1.49 

2.88 - 1.678 

6.91 -1.315 

10.90 - 1.634 

15.54 

1.00 

2.29 

4.93 - 1.082 

12.10 - 1.002 

19.35 - 1.000 

5.98 

0.75 

24.24 





0.50 







Similar frequency characteristics for er = 10 and the same 
different values of y considered before are presented in 
Fig. 2 and Table 1. The general arguments discussed for 
the case er = 1 hold. For y = 1.25, all the critical 
frequencies are decreased whilst the critical gains and the 
quality factor remain almost the same. For y = 1.00 the 


responses are hardly affected. For y < 1, all the critical 
frequencies decrease with decreasing y. 

Considering the phase characteristic, it is decreased 
(increased) for y = 1.25 (y < 1.00) and almost 
unaffected for y = 1.00. 
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Bode Diagram 



Table.5: Numerical data for frequency response curves when a = 10 for the filter. 


Y 

mi - 

W 2 (^2) 

a> p - (M p ) 

BW - ( BP ) 

Q - BP 

1.25 

0.537 - 8.815 

0.677 - 8.868 

0.618 -(11.844) 

0.677 - (0.140) 

4.414 

1.00 

— 

1.000 - (-3) 

0- (0) 

1.000 -(—) 

— 

0.75 

— 

0.145 -(-3) 

0-(0) 

0.145 -(—) 

— 

0.50 

— 

0.006 - (-3) 

0- (0) 

0.006 - (—) 

— 


Step response of the filter for a = 10 is shown in Fig. 5. But there is an apparent time lag (spread) for y = 1.25 > 

It is seen that the conventional Butterworth filter response 1, and time lead (squeeze) for y = 0.75 < 1. The lead 

is not affected by a when y = 1.00 (see the legendsy = (squeeze) is higher for y = 0.50. 

1.25 and y = 1.25 *). Further, the peak values of 
oscillations are not affected by changing erfrom 1 to 10. 
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Fig.7: Step responses of the transfer fiinction when a =10. 
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Table.6: Step response characteristics for a =10. 


Y 

T 

i r 

Tpi M pl 

T p 2 M p2 

Tp3 M p3 

T 

1 s 

1.25 

2.38 

7.23 - 1.678 

17.35 - 1.315 

27.37 - 1.634 

39.05 

1.00 

2.29 

4.93 - 1.082 

12.10 - 1.002 

19.35 - 1.000 

5.98 

0.75 

5.22 


— 


12.77 

0.50 

12.63 






VL CONCLUSION 

A FO Butterworth filter derived from the conventional 3 rd 
order Butterworth filter by the IO derivative to FO 
derivative transformation with two parameters y and o as 
in [17] is investigated in this paper. It is arrived the 
following conclusions. 

1. A variety of low pass filters can be obtained by 
changing the parameters yand o. 

2. For y > 1, the classical frequency response curve of 
the Butterworth filter disappears, and a peak occur 
in the frequency response. 

3. When y = 1. o does not effect on the 
characteristics. 

4. erdominantly affect the time characteristics, critical 
magnitudes of the time responses are not changed, 
but increase of a cuuses a slower (faster) response 
for y > 1 (y < 1). That is o effects like a time 
scaling operator. 

5. Effect of a on the phase characteristics is that 
increase of a decreases (increases the phase if y > 1 
(Y < 1). 

6. For y > 1 and as y -> 1.38 the filter can be used a 
high-Q narrow band pass filter as well. But in the 
limit case the filter becomes unstable. 

Having in mind these conclusions, fractional order low 
pass filters and/orband pass filers can be designed. 
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